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Abstract 

We discuss an extension of a map between between BPS states 
and free fermions. The extension involves states associated with a full 
two matrix problem which are constructed using a sequence of integral 
equations. A two parameter set of matrix model eigenstates is then 
related to states in SUGRA. Their wavefunctions are characterized 
by nontrivial dependence on the radial coordinate of AdS and of the 
Sphere respectively. A kernel defining a one to one map between these 
states is then constructed. 
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1 Introduction 



Studies [1-8] of giant gravitons in AdS Supergravity (and dual N=4 SYM 
theory) have lead to a simple (matrix model) picture for 1/2 BPS states. In 
particular a free fermion model [8,9,10] of harmonic oscillators was identified 
and shown to simulate fully the dynamics of 1/2 BPS states and their inter- 
actions. In [10] (referred to as LLM) a classical Ansatz for AdS (bubbling) 
configurations was constructed whose energy and flux were demonstrated to 
be in a one to one correspondence with those of a general fermionic droplet 
configuration. Further, relevant studies of this free fermion map have recently 
been carried out [11-33]. 

It is clear that it would be desirable to extend the map to more general 
states and go beyond the simple case of free fermions. This would require an 
investigation (and solution) of more complex two (or multi) matrix models, 
a formidable task. In the present work we present a step in this direction. 
We will attempt to extend the correspondence from the fermionic family of 
states (representing a single diagonal matrix) to a more general set associ- 
ated with states of a two matrix quantum mechanics. As was already seen in 
[8] which concerned itself with the case of 1/2 BPS states one can start with 
a system of two matrices, or a complex matrix and perform a truncation to 
a single hermitean matrix (in the manner analogous to a similar phenomena 
in the quantum Hall effect). The reduction was explained in [8] to be the 
Hilbert space equivalent of a holomorphic projection where the set of observ- 
ables are given by traces of the complex matrix Z only. The introduction of 
mixed traces, involving the second (conjugate) matrix immediately leads to 
a nontrivial dynamical problem whose eigenstates were never constructed. 

We will first adress this problem of constructing invariant eigenstates of 
the two matrix quantum system. For this we develop in some detail a hybrid 
formalism, treating one of the matrices fully in the standard collective field 
theory manner, while the other is treated in the coherent state representa- 
tion. This second matrix behaves then as an 'impurity'. The corresponding 
collective field theory of combined, mixed traces is then worked out and is 
shown to lead to a sequence of eigenvalue equations. These equations are 
seen to generalize an eigenvalue equation first found in and first solved 
for its eigenstates in describing angular degrees of fredom of the single 
matrix model. The sequence of eigenvalue equations can be solved for the 
present case of the oscillator potential. It provides a two parameter set of 
energy (dilatation operator) eigenvalues and a corresponding 2 dimensional 
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space of eigenfunctions. 

The central issue then becomes that of providing a correspondence be- 
tween the eigenstates of the matrix model and states and eigenvalues in 
Supergravity. Here we work in a linearized approximation specifying a class 
of fluctuations with matching quantum numbers. The wavefunctions, in the 
AdS X S background are nontrivial, being given by hyper geometric functions 
or corresponding special functions. Nevertheless, we describe a 1 — 1 map be- 
tween a two dimensional subset and the two dimensional set of wavefunctions 
given by the matrix model. This map involves a transformation introduced 
originaly in the context of the 2d Black hole and the corresponding matrix 
model jSnj. This transform, appropriately interpreted, then provides a one 
to one map between the gravity and matrix model wavefunctions. We em- 
phasize that being one to one this map is different from the well known 
holographic projection. It is expected that further studies of the map will be 
of relevance for reconstructing AdS quantum mechanics. 

The outline of the paper is as follows. In Sect. 2 we give a review of 
the simple matrix model and of the fermion map. In Sect. 3 we adress the 
two matrix problem describing its collective field formulation. We derive a 
sequence of eigenvalue equations and solve for eigenvalues and eigenfunctions. 
In Sect. 4 we consider the wavefunctions of the AdSxS SUGRA and specify 
the transform to the matrix model eigenstates. Several Appendices contain 
further details. 



2 Review 

We begin by reviewing and clarifying the existing map between the | BPS 
SUGRA configurations and the states of the harmonic oscillator matrix 
model. 

The matrix model degrees of freedom originate from a reduction of = 
4 Super Yang-Mills theory on R x S^. The hamiltonian is therefore the 
dilatation operator and the Higgs fields become quantum mechanical matrix 
coordinates $a (t), a = 1 . . . 6. For the study performed in the present paper 
one can concentrate on the dynamics of two matrices 

The commutator interaction did not play a role in the | BPS correspon- 
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dence and in what follows we will mainly concern ourselves with the simple 
quadratic harmonic oscillator model of two matrices 

H = -Ti {P^ + P| + $2 + $2^ 

The symmetries of this reduced theory are given by the U (1) charge 

J = Tr (Pi$2 - P2^i) 

and an SL{2,R) symmetry algebra (allternatively SU(2)). 
One has the complex matrices 

^($i+i$2) 
v2 

and the conjugates 



n = 
nt = 



1 


(Pi + iPs) = 


. d 






1 


(Pi - iP2) = 


d 


71 





Restriction to BPS configurations corresponds in the matrix model to con- 
sidering a subset of correlators involving only the chiral primary operators 
of the general form 

TrZ'=iTrZ'=2...TrZ'=". 

For the corresponding reduction in Hilbert space one proceeds as follows(see 
[8,9]). It is useful to introduce the operators 

A=]^{Z + iIi) 

and 

2 ^ ' 

In terms of these, the Hamiltonian and the U (1) charge read 

= Tr {A^A + B^B) 
J = Tr {A^A - B^B) 

4 



One now has a sequence of eigenstates given by 

Tr((At)") |0) E^J^n 
Tr((fit)") |o) E=-J = n 

Tr ((At)")Tr((5t)'") |0) E = n + m, J = n-m 

Restriction to | BPS configurations corresponds in the matrix model Hilbert 
space to a reduction to a subsector given by A osciUators. It is useful to 
diagonalize A, A'^ by using the unitary symmetry 

The measure in these variables shows that we can treat the Aj's as fermionic 
variables. The Hamiltonian for these fermionic oscillators is 

i 

The fermionic wavefunctions are 



iPf (Ai,A2,...,A„) = e-^'^'^'det 




x[^ . 
x^^ . 


• A^i^ 


\ 










J 



After dividing the wavefunction by the Vandermonde determinant we have 
that 

'^B;h,l2,...,lN (Al; A2, . . . , A„) = e ^''^''*''X«i,«2,.--,'iV (Al; A2, • • • , Atv) 

Where Xh.h,---,iN denotes the character of a representation of SU (N) that 
corresponds to a Young tableaux with /i boxes in the first row, I2 boxes in 
the second one etc. Of special interest is the sequence of states corresponding 
to representations that contains 1 row of I boxes 

ipB;h,i2,...,iN (Ai, A2, . . . , A„) = e-^^^'^'xz,o,...,o = e~^^^'^'X/A...,o (Ai, A2, . . . , Ajv) 

and another sequence that corresponds to a representation that contains 1 
column of / boxes 

^B;h,i2,...,iN (Ai, A2, . . . , A„) = e-^i^^^^xi,o,-,o = e-^*^'^»xi,i,...,i,o,...,o (Ai, A2, . . . , Ajv) ■ 
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In the fermionic picture [H7] the first set of states represents particles and 
the second holes. These were explained in [8,6,9] to corresponds to a giant 
gravitons in AdS and to a giant gravitons on the sphere respectively. In 
terms of the moments 

N 
j=0 

one obtains Schur polynomials representing these states 

Xl,0,...,0 (Al, A2, . . . , Xn) = Pi (01, 02, • • • , 0JV) 

Xi,i,-,i,o,...,o (Ai, A2, . . . , Xn) = (-)' Pi (-01, -02, • • • , -4>n) 

They are exact eigenstates of a cubic collective field theory representing the 
bosonized version of Id fermions. In terms of a two dimensional density field 
p {x, y, t) the hamiltonian is simply 

H = ^ J dx J dy{x^ + y^) p (x, y, t) . 

Together with the non-trivial symplectic form 

Lq = 2tc dx dyp (x) G {x — x) p (x) 



one has a topological 2 + 1 dimensional scalar field theory [28] which can be 
reduced to a 1 + 1 dimensional collective field theory describing the dynamics 
of the boundary (of the droplet) y± (x, t) by 

L = ^JdtJdx [y+d-'y+ - y-d-H- -{{vl- yl) + {y+ - V-))] 

One can parametrize the boundary in terms of radial coordinates, in which 
the Lagrangian becomes quadratic. This is a simple manifestation of the 
integrability of this theory. This goes as follows: 

Consider a closed curve r (s,t) in B? with parameter s which in our case 
describes the boundary of the fermi sea in the phase space. In general the 
equation of motion can be written in the form 

dtr X dsf = dgA (r) 
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with A (f) defining the model that we are studying. For the case of free 
fermions in an oscillator potential one has 



A(r) = -r". 
^ ' 2 

If we parametrize the curve as 

r (x, t) = XX + y± (x, t) y 

we recover 

dtr X d^f = dty± 

dty± = 'Id^yi + x') 

If instead one uses polar coordinates to parametrize the boundary 

^(0, = P (0, cos (0) X + p {(f), t) sin {(f)) y 
and in this case we have 

dtrx d^r = ^dtp"^ {(f), t) 
dtp' = d^p\ 

It is instructive to derive the above linear equation of motion from the non- 
linear one by using the field dependent coordinate tranasformation.lt is sim- 
ply given by 

X = p[(f) {x,t) ,t] cos {(f> {x,t)) 
y+ = p[<l>{x,t) ,t] sin {(f) {x,t)) 

These are then the action-angle coordinates for the dynamics of the 
boundary. 

In their work Lin, Lunin and Maldacena JUl have identified a nonlinear 
Ansatz for 10 d Sugra which exactly reduces to the above, bosonic hamilto- 
nian of Id fermions . To summarize the main features of the Ansatz, one has 
first the 10 dimensional metric 
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ydyVi 



h-'\dt + Vidx'f + K\dy' + dx'dx') + ye^d^l + ye 



-''dnl 



2y cosh G, 
- tanh G 



y{dy, - d,v;) 



and a corresponding Ansatz for the the gauge fields. 

The only unknown function z is shown to obey the Laplace equation: 



didiZ + ydy{^) -- 
which is solved boundary value problem 







z{xi,X2,y) 
Vi{xuX2,y) 



y"^ f z{x[,x'2,0)dx[dx'2 



TT [(x-x')2 + y2]2 



x')2 + ?/2]2 



Remarkably,the fiux and the energy of this general configuration were shown 
by LLM to take the form of the bosonized free fermion droplet 



N 



A = J dxi J dx2 [xl + xl) (t, xi, X2) + ^ 

J dXi J dx2 {xj + xl) {t, Xl, X2) + ^ 



It should be stressed that even though these expressions look two dimen- 
sional, effectively this is still only a 1 dimensional correspondence (it is de- 
scribed explicitely by the 1-t-l dimensional bosonic scalar field theory). In 
addition to the formulas for the fiux and the energy one also needs the sym- 
plectic form (which should coincide with the symplectic form established by 
Iso,Karabali and Sakita [^) for the 2d fermion droplet. Another, simple 
way to see the one dimensionality is by an analysis of linearized fiuctuations 
(we give this in Appendix A). One has 



S 



n>0 



dt 



_ -2 _|_ -2 



2 2 



in agreement with the well known quadratic action for chiral primaries in 
AdS: 



S^Yl / dx'Vgl^ [a--na^- - n (n - 4) a'^a^^] 

It is supersymmetry which requires that {dt — d^) cr = which for the + 1 
dimensional variables means that the "angular momentum" is equal to the 
energy. Choosing an opposite chirality for the fermions we would have had 
the condition {dt + d^) o" = which would flip the sign in the relation between 
energy and " angular momentum" . 



3 Matrix Model Eigenproblem 

We have seen in the discussion that the treatment of 1/2 BPS states corre- 
sponds to a reduction, namely to one matrix quantum mechanics given by 
the canonical set A and . It is our interest to extend this correspondence 
to a larger set of states. In the matrix model they will be states involving 
the two matrices {A and B) of a two matrix model. This can be stated as 
a two matrix problem, with two hermitian matrices M and N in a quadratic 
potential, i.e., with Hamiltonian 



H 



2 ^^dM dM 



1, 



d d 



) + -Tr{N^) (1) 



Using creation- annhilation operators for the matrix N^j in a coherent basis, 
the Hamiltonian takes the form considered in this article: 



18 8 1 8 
2 ^dMdM' 2 ^ ' ^ 8B' 



(2) 



We consider the action of this hamiltonian on functionals of invariant vari- 
ables (loops) 



^(A;,s = 0,1,2,...) 



where the s = 0, 1, 2, ...) are states with s "S impurities": 



9 



V'(fc,0) = Tr(e'^*0 

^{k,2) = r dk'Tr{Be"''^'Be'^''~'''^^) (3) 
Jo 

In terms of the eigenvalues Aj and the angular variables V of the matrix 
M = VAV^, we have 



^{k,l) = E^iV+BVhe'"^^ (4) 



^P{k,2) = -2zS„-(y+W),,(y+W) 



(A, - A.) 



Using the chain rule, we obtain for the matrix M kinetic energy operator on 
the wave functional: 



2 ^dMdM' 2 J ^ dMdM ' di){k,s) 

d^{k,s)dtlj{k',s'). 92 



^S,,,, y rfA; y rffc' rr(- 



dM ' dip{k,s)dip{k',s') 
As it is traditional ^39,, we introduce the notation: 



u{k, s) splits the loop ipi^k, s) and fl{k, s : k', s') joins the two loops 'ip{k, s) 
and ^/'(/c', s'). 

We will find it useful to introduce a density description, or x representa- 
tion: 
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/dk f 
— e"*''^?/'(fc, s), ilj{k,s)= dxe'''''ilj{x,s). 

Any function of k (or x) transforms accordingly. Namely: 

-ikx, 



uj{x,s) = / — e uj{k,s) 

ZTT 



nix, s;y,s') = J^l^e-^'^e-^>''ynik,s;k',s') 
For conjugates, we have 



dip{x,s) J dijj{k,s)' dijj{k,s) J 2tt dijj{x,s) 

In the density description, the kinetic operator then becomes: 

-is.,,. / / . : 



3.1 Spectrum and fluctuations in the zero impurity 
sector 

Consider first the analysis for the spectrum of the zero impurity problem. 
This sector corresponds to the Quantum Mechanics of a single hermitean 
matrix, and it has by now a standard solution jHlI , 110] , IS] , which is briefly 
reviewed in Appendix B. In this case, one has the standard cubic Hamiltonian 

H!ff = dxdM^)ijix,0)dM^)+N'(^ J cixy^3(a:,0)+^(x,0)(y-^) 

(7) 

giving the well known Wigner distribution background in the limit as — > cxd 



TTip{x, 0) = 7r0o = \/2[i — x'^ = V2 — x'^. 
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For the small fluctuation spectrum, one shifts the background 
to find the quadratic operator 

^2 = ^ y dxTK^oP^ix) + dxTlMdxVf 

The way to diagonalize is by now well known: one changes to the classical 
"time of flight" q 



— = 7r0o; x{q) = -V2cos{q); 7r0o = v^sin(g); < g < vr 
aq 

One obtains the equation for a 2d massless boson: 

^2 = 1 J dqP\q) + IJ dq{d,^f (8) 

In addition one needs to impose Dirichelet boundary conditions at the 
classical turning points, for a consistent time evolution of the constraint 
fl52|l ■ Therefore the spectrum in the zero impurity sector is 

Wj = j ; = sin(jg) (9) 

The following comment is in order: the harmonic oscillator potential is 
special, in that the effective hamiltonian ((7j) can be equivalently written 
as (for discussions on the relationship between the two re-writings in the 
context of supersymmetric or stochastic stabilizations, see for instance 



H^eff = dxdM^)ij{x,0)dMx) + ^ J dx^{x,0){ j dy^i^^-x) 
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It is then seen that the Wigner distribution background also satisfies the 
well known BIPZ ^H] equation 

dz f''^^ ) _ ^ /-|^Q\ 
[x — z) 
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Shifting about the background as above, we obtain for the quadratic hamil- 
tonian 



This non-local hamiltonian can be easily shown to be equivalent to (jH)). Let 
us examine this in slightly more detail: by changing to the classical time of 
flight q, we obtain 

The above non local integral operator plays a prominent role in what follows 
and is discussed in Appendix C. Let us denote it by 

and by abuse of language (it does not satisfy a Leibnitz rule) refer to it as 
the "absolute derivative", for ease of notation. We note that {—i\dq\y = 
and that the appropriate eigenfunctions of this operator are 0„ = sin(ng) 
with eigenvalue n as shown in Appendix C. Therefore the eigenfunctions (jHI) 
are also the solutions of 

{idt + t\d,Mq) = 

3.2 Quadratic hamiltonian for states with impurities 

We return now to the (pre-hermitean) kinetic energy operator ^ (or 
We note that 

< ^p{x, s) >=< ^lj{k, s) >= ; s = l,2,3,... 

This observation implies that for the multi-impurity spectrum it is sufficient 
to consider the zero impurity sector Jacobian already discussed jUj, i.e., 

J^-z — : rJ 2 = -- — — -In J 



dilj{x,0) di/j{x,0) di/j{x,0) 2d^lj{x,0) 

d d 

s = 1 2 3 ... 

dipix, s) dip{x, s) ' ' ' ' 



13 



where, to leading order in 



Let us now identify the terms in © which determine the quadratic operator 
in the multi-impurity sector. 

We look for terms of the form ipi^y s)d/dilj{x, s), s > when ipi^^ 0) — (poi^)- 
Contributions of this form contained in the first term of (jH)) result from 
splittings of the loop ip{x, s) into a zero impurity loop and another with s 
impurities. We will denote this amplitude by lj{x,s). 

Contributions contained in the second term of © are obtained as a result of 
the similarity transformation described above, when we replace d/dip{x, 0) — >■ 
— {l/2)d/d'ip{x,0)\n J . We therefore obtain: 



If d Iff dlnJd 

if| = — - / dxLD{x, s)——- + - dx dyfl{x, : y,s) 



2 J ' d^jj{x, s) 2 J J ' • ' ^^p{x, 0) d^{j{y, s) 

(12) 

In a problem involving joining and splitting of loop states, the issue of 
closure of loop space is an important one. The first term in ()12|1 always 
closes. This is because 

u{k,s) = -2 dk'k'^{k\s)i){k-k' 
Jo 

This result is a straightforward application of a result established in |1S] . In 
the X representation. 



dk 

^{z,s) = / —e-''''u:{k,s) 



il){z,s) I dx 



%[x) 



ip{x, s 



(poiz) / dx- ^^ + / dx(f)o{x)d^{- 



[X — z 



ij{z,s) 



[Z — X) 



Substituting this expression into (fT^ we obtain 
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^ ' J {x — zY dijj{x, s) 

r M^)ij{x,s) d 

'^''J'^' (x-z) '^'di:{x,s) ^^'^^ 

1 /•,/•, , dlnJ d 

+ 7: dx dyU[x,0 : y, s) 



dip{x, 0) dipiy, s) 

In general, for an arbitrary potential, the last term in (fT^ involving fi(a;, : 
s) will not close. We will argue in the following that for the harmonic 
oscillator potential this term closes, by considering explicitely s = 1, 2, 3 and 
then arguing for the general case. 

3.3 The one impurity sector 

It is straighforward to show that in this case 

from which it follows 

n{x, 0; y, 1) = d^dy{ilj{x, l)S{x - y)). 
The term involving VL[x, : y, 1) in (|13|) becomes 



\ j '^^ j dy^ix, 0:y,l) 



din J d 

dil>{x, 0) dil){y, 1) 



'(9'?/'(x,0) ' dil){x^l) 
dx I dz f'^^ \ iIj{x, l)dx- ^ 



{x — z) ' dip{x,l)^ 

where we have used flll|) . We observe that this term cancels exactly a similar 
term in (fT^ , and we obtain the final form for the quadratic hamiltonian in 
the 1 impurity sector: 

Hr' = [dx[ d /o{zMx,l)-i.iz,l)Mx) d 

J J [x- zf di)[x, 1) 
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The rescaling (j55p leaves the above hamiltonian invariant, or equivalently 
the above hamiltonian is of order 1 (N^) in A^, as was the case in the zero 
impurity sector. Writing the operator as 

f d 

dx i dyilj{x,l)K{x,y)- 



we obtain 

Acting on a wave functional 



we obtain the Marchesini-Onofri kernel . |35j . |49j . [HO] 

Using (fTUI) . the first term yields —f{x), and by changing to the time of flight 
coordinates the kernel can be written as 

or, for the spectrum equation 

(-1 - ^|9,|)(7r0o(g)/(g)) = winMq)fiQ))- 

As described in Appendix C the spectrum and eigenfunctions of this operator 
are 

1 sinfrig) 

Wn = n-1 ; 0„ = ■ n = l,2,... 

y/Z sm(gj 

For the the harmonic oscillator potential, these are the well known Tcheby- 
chev polynomials of the second kind. Adding the contribution from the 
Tr{Bd/dB) term of the Hamiltonian we obtain 

Wn = n ; </>r^ = ^^^ ; n = l,2,... (15) 
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3.4 The two impurities sector 

For two impurities, we have 



-2koi:i,jiV+BV)ijiV+BV)j^ 



— ik 



(A. - A,) 



+ e 



ikoXi 



(A. - A,)2 



and 



n{x,0:y,2) 



-2zd.dy^.,{V+BVUV+BV),Mx - yVj^ 
-2zd,j:,,{V+BVUV^BV),,6{x - A,)fcg^|^ 



The fl{x, : y,2) term in ()13p takes the form 



—i dx dy 



^,,{y+BVUV+BV),,5{x - y)^-}^L^d^^^hl^d^, ^ 



- Y.,^,{y^BV),^{y^BV),,b{x-\ 

= -2i J dxJ:^J{V^BV)ij{V^BV) 



(A, -A,) ^d^ix,0) "d^iy,2) 

S{y-\,)-5{y -\,)^ d\nJ d 

(A. - A,)2 
6{x — Xi 



^^p{x,0)^^p{y,2) 



(X - A,; 

2z / 



{x — z) dip{x, 2) 



dz 



(A. - z) 



dz- 



d 



{X,-z)\di^{y,2) 



For the harmonic oscillator potential, we can use the result fjlO|) . so that 



91n J 



= -2i J dxi:i,jiV+BV),j{V+BV) "^(^"^^^ 



+2i / dyT.i^j{V+BV)i,{V+BV),, 



^^'x — \j) J "'~ {x — z)^'^ dilj{x,2) 
6{y - A,) d 



(y-X,) d^P{y,2) 



(16) 
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But from (g)), 



/flh 
^e-^'-^Pik 2) 
Ztt 

/flh p^^^ 
—e-'^^i:,,{V+BV\,{V+BV) 



5{x — \j) 



{x - Xi) 

This allows us to express entirely in terms of the density ip{x, 2) as 



1 I dx I dy^l{x,0 : y,2)- 



^'^|J{x,0)^'^p{y,2) 



dx I dzp^i){x,2)d^—^-— - [ dxij{x,2)- ^ 



{x — z) ' dip{x,2) J ' dil){x,2) 

As was the case for the one impurity sector, the first term above cancels the 
similar term in ()13|). and we obtain for the quadratic hamiltonian in the 2 
impurity sector: 



rrs=2 [, f, 0o(^)^(x,2)-^(2:,2)0o(x) d f d 

(17) 

This is a shifted Marchesini-Onofri operator. It can be recast in the form: 



(-2 - t\d,\){nMq)fiq)) = w{iTMq)fiq))- 

The spectrum and eigenf unctions of this operator are 

,_9 sinfrio) 
Wn = n-2 ; 0r = • )\ 5 n=l,2,... 

V2sm(gj 

Adding the contribution from the Tr[Bd/dB) term of the Hamiltonian we 
obtain 

Wn = n ; 0r'= ^ . :\ ; n=l,2,... (18) 

V2 sm(g) 
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3.5 Multi-impurity spectrum 



The pattern that emerges from the above discussion is clear: for s impurities 
and the harmonic oscillator potential, one obtains a shifted Marchesini-Onofri 
operator with spectrum and eigenfunctions 

sin(ng) 

Wn^n-s ; 0^ = ——— ; n=l,2, ... 

V 2 sm(g) 

When the contribution from the Tr{Bd/dB) is added, we have for the full 
hamiltonian 

wn^n ; 0:-^^^ ; n = l,2,... (19) 

To provide further evidence of this pattern, the 3 impurity case is treated 
explicitly in Appendix D. We also checked that by introducing multi local 
densities and then projecting to the 2 and 3 impurity states discussed here, 
we obtain the spectrum described above. 

To summarize, as the U(l) charge operator J it is represented by 

and consequently j = n — 2s. Together with the energy eigenvalues w = n, 
these specify a two parameter family of states and a two dimensional complete 
set of eigenfunctions. 



4 SUGRA Map 

In this section we would hke to identify the states of Sugra fluctuations and 

establish a one to one map with the eigenstates of the matrix problem found 
in the previous section. With the two matrices we hope to explore the extra 
coordinate which will be related to the radial coordinate of AdS and S. Since 
the other angular coordinates are ignored, it is sufficient to concentrate on 
the small fluctuation equations asociated with AdS^ x 5*3 (the analysis for 
AdS^ X 5*5 reaches an identical conclusion). Wc have obtained in the matrix 
model solution a two parameter sequence of states with the eigenvalues J = j 
and w = j + 2n. It is easy to find a corresponding sequence of states, 
which have the same eigenvalues. Actually there are two sequences ,one with 
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nontrivial functional dependence in the radial variable of ^4^5' and the other 
in S. This situation is familiar from giant gravitons. 

It will be clear that while the integer valued eigenvalues easily agree 
(between the matrix model and supergravity) , the comparison of their wave- 
functions is much less trivial and also much more interesting. In Sugra the 
wavefunctions are given as nontrivial special functions, while in the solution 
of the matrix eigenvalue problem they take the form of ordinary plane waves 
. The later obviously happens after the change from eigenvalue coordinate to 
the "time of flight" coordinate. We will establish a relationship between the 
two pictures in terms of a kernel describing a (canonical) change of variables. 



4.1 The LLM kernel 

It is useful first to work out the form of the kernel for the case of 1/2 BPS 
states given by the LLM map. For this one has to consider the LLM construc- 
tion and perform the small fluctuation analysis . We do this in Appendix A 
where we also give the details of a transformation to the Lorentz-De Bonder 
gauge. Furthermore, from now on the time of flight q will be denoted by r. 

Let us concentrate on the fluctuations associated with the metric g^^. In 
the gauge of LLM the perturbation Sqqq reads 



'''' y I - ^UAdS {1 + 2uAdsf 



sin^^— / dT ^^^^^ oVa.e^^- (20) 

27rJo [l + a2-2acos(T-</.)]'^ ' 



cos^ 

a 



cosh p 

The relevant gauge transformation can be written in integral form as 



59^ 1^'"'^ Ta,e-^ 

cosh p — cos^ ^-^ 



3 



2 .2. , o (21) 



= — tant^— / dr- — — > 

2nl-a^ Jo 1 + a2 - 2acos (r - 0) ^ 
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Performing the gauge transformation we have 

%n = sin^ 9- / dr \ " ^ ^ a,e^^- (22) 

27r Jo [1 + - 2a cos (r - 0)J 

In this form we see the relation 

, , , , 1 /"^^ 1 - - + cos (r - 0) • • , , 

2 1 J I a, (t, pA,e) = - dr \/' Yl ^3^' (23) 

27r Jo [1 + - 2a cos (r - 0)J ~^ 

After performing the field dependent gauge transformation in order to 
recognize the primary field coming from the metric and the three form one 
has the relation 

\j\ {t, p, 0, 9) = — e»^* drK'^'^'' {p, 0, 9\t) e'^^ (24) 
where the kernel is given by 



[l + a2 -2acos (r - 0)]^ 



cos^ 

a — 



(25) 



cosh p 

At this point we notice that a < 1 at points where the measure of AdSs x 5"^ 
is non-zero. At this places we introduce a cut-off L limiting the angular 
momentum j .We then have the kernel 



1 r 

\j\ aj {t, p, <P, 9) = -e'^* drKt'^'' {p, <f>, 9\t) e'^\ \j\ < L 
= ^e'^' drK^'^'' {p, 0, 9\t) e'^\ \j\ < L, \j\ > L. 



(26) 



The kernel with the cutoff is given by 
i^f^^^(p,0,^|r)=X™(p,0,^|r) + 
— cos [L (r — (f))] + a cos [{L — 1) (r — 0)] ^ 

1 -I- a^ — 2a cos (r — 0) 
aL cos [{L + 2) (r - 0)] - [1 + L + 2La'^] cos [{L + 1) (r - 0)] ^ (27) 

[l + a2 -2acos(r-0)]^ ^ 
-a^ (L + 1) cos [(L - 1) (r - 0)] + a [2 + 2L + La^] cos [L (r - 0)] ^ 
[l + a2 -2acos (r - 0)]^ " 
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We see that we have a strong convergence 



L—*oo 



4.2 Correspondence with the 2d Black Hole 



(28) 



To proceed with the construction of the kernel in our more general two di- 
mensional case it is also useful to take note of a correspondence with an 
equivalent problem that was considered in the case of a 2d black hole. We 
show in what follows that there is a simple connection between "off-shell" 
black hole wavefunctions and on-shell AdS wavefunctions that we have iden- 
tified. 

The wavefunctions that we consider correspond to highest weight states 
on SO (4) but with a nontrivial dependence on the radial coordinate of AdS 



/ = cos' (0) e'^^i) {t, a) 
We have the following eigenequatin for ip 

- cos^ (a) dill) + cos^ (a) dl + cot (a) d^^p = l{l-2)^p 



cos^ (cr) 



COS (cr) sin (cr) 
with the integration measure 

^^g^^dtda = tan (cr) da. 



dm 



A change to a new function 



R 



cos (cr) 



with the new measure 
leads to the equation 

-^R+^^ 

cos cr 



dm — - sin (2cr) dtda 



dl + 



1 



cos (cr) sin (cr) 



cos (cr) R 



/ (/ ~ 2) 



cos^ cr 



R 



or 
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dlR + 2 cot {2a) d^R 



/ (/ - 2) + 1 



cos^ cr 



R~uj'R + R 



This can be compared with the 2d black hole equation [Sd] defined as a coset 
SL (2, R) /U (1). For the case of the Lorentzian black hole they are specified 
by the eigenvalue equation 



1 



d^^T^ + dlT^ + coth (r) drT^ 



X' T. 



4sinh(0 " 

— ^Z/2yA ^ g2j.X ^ ^^^^ g^j.^ 

smh (2) 



A 



-- - AM T 



and the inner product is defined through the integration measure 

/■oo 

Jo 

We see that the two problems are related, through the following transfor- 
mation transformations 



a 



- (r + Tc) 
2^ ' 



In ref. jHS] a transformation was constructed relating the wavefunctions 
in the black hole case to those of a c=l matrix model. The transformation 
reads 



dtn 



da5 



sinh 



sinh 



2tn 



cosh {2a) 



e-'^^ cos (4Acr) . 



and it involves a nontrivial kernel which specifies a canonical transformation 
from one problem to another. In the present case we will follow the construc- 
tion of and construct an analogous kernel which will relate AdS (and S) 
wavefunctions to those of the matrix eigenvalue problem. 
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4.3 The AdS Kernel 



We first give the main formulas defining the AdS kernel. The wavef unctions 
obey the equations 



(29) 



Visa,- „ (t, p, (f),e) = - \j\ (IjI + 2) a,- „ {t, p, 0, 9) 

^AdS3^j,n {t, P, 0, 6) = \j\ - 2) CTj- „ {t, p, 0, 6) 

d 

- ^'o^^hn {t, P, 0, 6*) = jcrj,n {t, p, 0, 0) 

and have an explicit solution in terms of hypergeometric functions 



aj,n {t, p, 0, e) = em^'^^>* cosl^l OeWl'^'l' cosh-(l^l+2") pF {l - j - n, -n; 1; - sinh' 
'j\+2n 



We now use the integral representation 



dz 

cosh p 



coshpe IjI"'^ (coshp + 2;sinhp) 



X 



+ sinh p ] (cosh p + z sinh p) e ^ i^i 
^ / 



1 



c 



i2TTZ 



- Iil+2n 






e i-'i V 



where 



w 



V — 



cosh p (cosh p + -2 sinh p) 
/ cosh p 



V ^ 



+ sinh p 1 (cosh p + z sinh p) 



bl+2n<L 
to derive the kernel defined through 

-I /'27r /»2-7r 

= em-^>*^ d^y^ drii'^(p,0,^|a,r)e^m»^")-+-], 
i^L (P, 0, e\a, t) = I ^ (Fl {w\t) Gl {vW) - 2Fl {w\t) Gl (via] 
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(30) 



(31) 



(32) 



(33) 



(34) 



Explicitely, the functions involved in the definition of the kernel can be 
worked out after an introduction of a cut off L for convergence. They take 
the slightly long forms: 

^ , , , 1 — Aw"^ — w^ + Aw^ cos (r — 0) 

Fl (w|t) — ^ — 2 ^ 

+ — 2w cos (r — (f))] 

— cos [L (r — (J))] + w cos [{L — 1) (r — 0)] 



1 + — 2w cos (r — 0) 
«;L cos [(L + 2) (r - 0)] - [1 + L + 2Lm;2] cos [{L + 1) (r - 0)] 

[1 + — 2w cos (r — 0)]^ 
-w^ {L + 1) cos [{L - 1) (r - 0)] + tt; [2 + 2L + Lw^] cos [L (r - 0)] ^ 



[1 + — 2-u; cos (r — 0)]^ 



1 



+ f 2 — 2t> cos {a — 20) 
- cos [L {a - 20)] + V cos [(L - 1) (a - 20)] 



(35) 
(36) 



Fi(«;|T) 



l + v^-2v cos (a - 20) 
1 + w'^ — 2w cos (r — < 



— cos [L (r — 0)] + w cos [(L — 1) (r — 0)] ^ 
1 + — 2^ cos (r — 0) 



Gl (via) = — ^ ^+ 

^ ' ^ [l + T;2-2^;cos(r-20)]' 

vL cos [(L + 2) (r - 20)] - [1 + L + 2Lv'^] cos [(L + 1) (r - 20)] 
[l + 'i;2 _2vcos (r-20)]^ 



(37) 



-V 



.L 



-v^ {L + 1) COS [(L - 1) (r - 20)] + v [2 + 2L + Lv^] cos [L (r - 20)] ^ 
[1 + 'u2-2^;cos(t-20)]^ ^ 



(38) 



4.4 The Sphere Kernel 



We now consider the second sequence of wavefunctions,which are character- 
ized by a nontrivial dependence on the radial coordinate of the sphere. The 
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wave equations read 



Vlsaj,n {t, p, <f>,e)^-{\j\+ 2n) + 2n + 2) a,- „ {t, p, 0, 9) 
^Mss(^j,n (t, P, 0, ^) = (b1 + 2n) i\j\+2n- 2) (7,-„ {t, p, 0, 9) 

d 

- ^^^hri {t, p, (j), 9) = jaj^n {t, p, (j), 9) 
In the coordinate system where the metric is 



ds^ = - cosh^ pdt"^ + dp"^ + sinh^ pdi/j"^ + d9'^ + cos^ 9d(j)'^ + sin^ 6'#^ 
the normahzable solutions are given by 

aj. „ = e^^'^^-"*e*^'^cosh-l^l-2"pcosl^l 9F (l + |j| + n, -n; l]sm^ 9) 
i^j,n =\j\ + 2n 



(40) 



We use the integral form of the wavefunctions 

dz ( cos 9 — z sin 



ii^j.nt ij4> 



i27rz 



cosh p 



cos^ + siii^ 

z 

COS ^ — 2; sin ^ 



/ cos 6* — zsin^ ,■ 



c 



i27iz 

dz 
i2'Kz 



cosh p 



|i|+2n 



cos6' + 



sin 9 



COS ^ — 2; sin ^ 



e 



\j\+2n 



e W 



where C is the unit circle on the complex plane and we defined 

cos 9 — z sin 9 



w 



V — 



coshp 
cos^ + ^ 

z 



(41) 



(42) 



cos ^ — z sin ^ 

Introducing a cut off on the angular momentum as we had for the LLM 



case 



+ 2n<L 



(43) 
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we rewrite the wavefunction as 

c2n r^TT 



j 1 f f J 

+ 2n) (J,- „ = em--*_ J da J (ItKl (p, 0, e\a, r) e'^ 
r dz 

KL{pA.O\a,T)= i --FL{w\T)GL{v\a). 
Jq iZnz 

(44) 

The functions {w\t) and (f |cr) specifying the kernel in this case are 
found to be given by 

^ , , , 1 — 4w^ — + 4:W^ cos (r — d>) 

Fl {w\t) — "2 — - + 

[1 + w"^ — 2w COS (r — 0)] 
— cos [L (r — (f))] + w cos [{L — 1) (r — 4>)] 
1 + w"^ — 2w cos (r — 0) 
cos [{L + 2) (r - 0)] - [1 + L + 2Lw2] p + 1) (r - 0)] 
[1 + — 2w cos (r — 0)]^ 

L + Lw^] cos [L (r - 0)] 

(45) 



-w^ (L + 1) cos [(L - 1) (r - (/))] + w [2 + 2L + Lw'^] cos [L (r - 0)] 



[1 + — 2^ cos (r — 0)] 
and 



1 - 1;' 

Gl {v\a) = 



l + i;2-2vcos ((T-2(/)) 

- cos [L {a - 2q')\ + ('COS [(I - 1) (a - 2c>)]^l 



(46) 



l + t;2-2vcos ((7-20) 

Let us now make the following comment regarding the cutoff that we have 
used. Since it imposes an upper limit on angular momenta it clearly plays 
a role of the 'exclusion principle'. Its removal seems to lead to singularities 
both in the sphere and the AdS case. One should remember then that this 
analysis is done at the linearized level, so there is no essential difference 
between the two cases. We can also show that if we restrict our attention 
to 1/2 BPS wavefunctions (which would correspond to n = 0), the above 
kernel reduces to the kernel that we have found from the LLM construction. 
We notice that the function Fl {w\t) is analytic in the unit circle C of the 
z— plane for every value of the remaining variables. 

2 /'27r /»27r 

|J I (^jfi = e^'*^ dr daKL (p, 0, e\a, r) e'^^ (47) 
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Performing the integral over cr gives 



/ daGLiv\a) = l (48) 
Jo 



which estabhshes the result 

f27r 



f C dz 

/ daKL{pA,e\a,T)= i Fz.(w;|r) = Fi(«;|r)|^^o = i^™(p,0,^|r) 



(49) 



5 Conclusion 

We have in the present work considered the simple a complex two matrix 
model with a purpose of developing further its correspondence with AdS 
eigenstates. We develop a (hybrid) formalism to construct a two dimensional 
sequence of invariant matrix model eigenstates. Here one of the (matrix) de- 
grees of freedom is treated in a density representation (in a manner analogous 
to the one matrix collective field theory), while the other is represented in 
the coherent state picture. This leads to a sequence of (integral) equations 
which we then solve for the case of the oscillator potential. The two dimen- 
sional set of eigenstates extends the one dimensional space representing the 
eigenstates of free fermions. As such this extension allows a nontrivial probe 
of one further extra dimension . This as we argue can be mapped into either 
the radial coordinate of AdS or the radial coordinate of the sphere. 

The mapping between states of the matrix model and the wavefunctions 
of SUGRA is one to one. As such it differs from the holographic map where 
one of the dimensions is projected out. In the present case the map can be 
described by a (two dimensional) kernel in paralel with similar maps in the 
case of 2d noncritical string theory. We also note that leg factors of this kind 
were found in the pp-wave map of |54j. 

When applied to a one dimensional subspace of 1/2 BPS wavefunctions 
our kernel reduces to the (linearized) map of LLM. In the construction of 
the extended map one seemingly requires a cutoff providing an interesting 
implementation of the 'exclusion principle'. The understanding of this cutoff 
is clearly of further interest. 

It should be commented that much like in the 1/2 BPS case of free 
fermions the model considered is that of simple decoupled harmonic oscil- 
lators. Yang-Mills type interactions present in the full theory might be of 
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relevance but are not included in our study. For the case of 1/2 BPS cor- 
relators there are theorems regarding the absence of coupling constant cor- 
rections. It can be hoped that this will persist for the present set of states. 
Certainly, the effect of coupling constant correction deserves to be investi- 
gated (e.g, |55j|)- It is also of interest to extend the present map to a still 
larger set of eigenstates. 
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7 Appendices 



7.1 Appendix A: Expanding the LLM solution in fluc- 
tuations 

In this section we would like to expand the circular droplet solution in "off- 
shell" fluctuations of the matrix model and see the equations of motion these 
fluctuations satisfy from the bosonic equations of motion of gravity. This anal- 
ysis was also performed independently in a recent paper j2'6j . 

The general iBPS LLM solut ion for the metric is determined by the 
function 

u{xi,X2,y) = — / dx'^u {xi,X2,0) 



TC 



X — X ] + 



with u {xi,X2,0) being the phase space distribution of the fermions in the 
matrix model picture. Parametrizing the boundary of the fermi surface using 
the polar coordinates representation 



xl (0, t) + xl (0, t) = R\^s + X]^" ^^'^^ + ^^'^^ 

n>0 

the phase space density becomes 

u (r, 0, 0, t) = -9 ^^R\ds + YlPn (t) sin (n0) + ng„ (t) cos (n0) - j + ^ 

and after approximating at first order in pertubations pn, Qn the distribution 
becomes. 



u (r,0,O,t) 



2 

AdS 



r) 

E 



,n>0 "^^AdS 

The field that is produced is given then by 

u (r, 0, y, t) = UAds (r, 0, y,t)+u (r, 0, y, t) 



Vnit) . . ,s , quit) 

sm [n(h) + — cos [nc, 



'^R-AdS 



u (r, 0, y, t) = UAds {r, 0, V^t) - — 



2 . I]„>oPri (t) sin (?20) + ng„ (t) cos (n 



RAds + r' + y'- 2rR\,, cos - 
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The above integral can be computed from the more general 



f27r ^im4 



f 

Jo (a — 2 cos (0)) 



- (h dz- 



c 



iz"^ -az + l) 



i (/) dz- 



z - z+f {z - Z-f 



27r 



d 



z± 



dz [z - z+) 
a ± ^/a? - 4 



= 27r- 



{z- - z+y 



n + 



z+ + z. 

Z^ — Z- 



Where the contour C is the unit circle on the complex plane of integration 
and we only picked the contribution from Z- which is the root that is inside 
the circle for a > 1. 
After setting 



y 

r 



R^Ads ^^"^^ P sin d 

?2 



^Ads cosh p cos 9 



the result is given by 



u {p, 4>, 9, t) = UAdS {p, 4>, 9, t) - 



sinh psin^^ 



^AdS 



E 

n>0 

where 



cos^ 



cosh p 
UAdS (p, 0, 0, t) 



n + 



cosh p + cos^ 9 
cosh^p — cos^ 9 

1 sinh^ p — sin^ 9 

2 sinh^ p + sin^ 9 



(cosh^ p — cos^ 9) 
\pn {t) sin {n(f)) + nqn {t) cos {n(f))] 
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The perturbation of the metric on is given by 

^AdS 

(cosh^ p sin^ 9 + sinh^ p cos^ 6) — i^M^=ud9'^ 

sinhpsin6' V 1 - 41*^^5 

4 cosh p sinh p sin^ OV^'^^ ~ 

4 sinh p sin 8 sinh p sin ^i^AdsK^'^^^ 2 cosh V cos^ ^^Ads ~ 
\4 H ^ u H „ u 



a/I - yr^^M^^ sinh p sin ^ y/l - 4.v?j^ 



AdS 



— 2 smh psni9\ I „ '^^'^ 7yUdVt\. 

l-2uAds{l + 2uAdsf 

At this point we would hke to show that the degrees of freedom g„,Pn 
turn on the chiral primary fields of IIB SUGRA on AdS^ x S^. After 
performing the field dependent coordinate transformation 

1 sin ^ cos ^ / COS^\"', / \ . / ,\ /\ / ,\n 

TTBi To T7i TT \Pn {t) sm + nqn [t) cos (n0) 

2^1d5 cosh" p - cos2 9 \coshpJ 



9^9- 



n>0 

1 cos^^^tanhp f cos^\", / x • / ,\ /\ / ,\t 

^ ^ ^ + ^ coshV-cos^^ g J ^" ^"'^^ + ^"'^^^ 

^y^(— T— ) bn (^) cos (n0) - g„ (t) sin (n0)] 
Ads „-^n Vcosnp/ 



^'^'S' n>0 



the 99 and 6*3 components of the first order perturbed metric are scaled by 



[n 



+ 1) ( cosh^ ) ^""^^ ""^^ ^'^^^^ ' 



d4 

After this observation we may identify the chiral primary fields 



8i^Ad5 ^ Vcoshp 



The correctly normalized action for the chiral primaries as given by Seiberg 
et al. reads 

S = Y1 ^2 [ dx'V^ [a-"na+" - n (n - 4) a-a^"] 

~~ in+1) jAdS^ 
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which after performing the spatial integral on AdS^ gives 

-2 I -2 2 2 2 

—Pn + (ln-n Qn-Pn 

and for each n we have a four dimensional phase space. Supersymmetry 
requires that {dt — d^) cr = which for our + 1 dimensional variables means 
that the "angular momentum" is equal to the energy. Choosing an opposite 
chirality for the fermions we would have had the condition {dt + d^) a = 
which would flip the sign in the relation between energy and "angular 
momentum" . 

7.2 Appendix B: Hermiticity and the zero impurity 
sector 

The zero impurity sector is the usual single matrix problem for Mij. We are 
interested in fluctuations about this single matrix background. As is now well 
known [SH] , 110] , IH] , this background is only exhibited as the stationary point 
of an explicitly hermitean effective potential. We recall the construction of 
this effective hamiltonian|39j. 

In order to take into account the non trivial Jacobian J involved in the 
change from the original variables to loop variables, one needs to implement 
the similarity transformation {i is a generic loop variable) 

di J^diJ~^ = di — -^Sjln J 
The Jacobian satisfies jSHl 

^ijdj In J = uji — djVlji 

The terms of the kinetic energy operator that are sufficient to generate the 
background and fluctuations are then [SIl , , |S3 

1 1 

--di^i^dj + -uJi^^ujj (50) 
In the zero impurity sector, 

cu(fc,0) = -k dk'i){k',Q)i){k-k',Q) 
Jo 

n(A;,0;fc',0) = -kk'tp{k + k' ,0) (51) 
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The X representation of ipik, 0) is the usual density of eigenvalues: 

0) = - Aj), 



and 



O(x,0;?/, 0) = d.j,dy{i){x,Q)5{x-y)) 

^(z,0)' 



u;(x,0) = -2d^{i){x,Q) j dz"!^^^^ 

From (j^Uj) we then obtain the form of the effective hamiltonian which is 
sufficient for the discussion of background generation and fluctuations: 

+ y"dx7/>(x,o)(y-/i)) 

where the Lagrange multiplier enforces the contraint 

j dx^{x,0) = N. (52) 



Since 

d^dyn~'^{x,0;y,0) 



S{x - y) 
^{x,Q) 



and 



the effective hamiltonian becomes: 



1 /■ , . d . d 

dxOx 



(54) 



2 J "d^{x,0)'^' ' ' dilj{x,0) 
To exhibit explicitely the N dependence, we rescale 
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X VNx 
^p{x,0) ^VNtlj{x,0) 

and obtain 

H!ff = dxdM^)^Pix,0)d,n{x)+N^(^ j dx^^3^a;,0)+V'(x,0)(y-/i) 

(56) 

which is equation ((7j) in the main text. 

7.3 Appendix C: Marchesini-Onofri Kernel 

We consider the problem of finding the spectrum of the operator 



(87) 



We start with the second term and consider the following integral, in "time 
of fiight" coordinates: 

— 7r0o(g)^^ ^>0 

_^ TT x(go) - x{q) 

Note that the range of the integral extends over a full period 2L = 2tt of 
the classical motion. Therefore, the integral above can be calculated by the 
residue theorem, by choosing a vertical path from —n + ioo to — vr, then along 
the real axis from — vr to vr, and then along a vertical path from vr to vr + ioo, 
"closing" at +ioo. The contribution fom the vertical paths cancel, due to 
the periodicity of the classical motion. The origin of the "time of fiight" can 
always be chosen so that the only poles on the real axis occur at g = go 
and q = —qo, corresponding to an even (in q) "displacement" x{q) and odd 
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"velocity" ircpQ^q). We alwyas choose a principal value prescription for poles 
on the real axis (half of the residue). If there are no other poles, as is the 
case in general for stabilised potentials, we obtain the result: 



r dq , , , e^"« ^. r dq , , , sin(ng) ^ . , , 

/ — TT^oig)^^ ^ = 2t / — 7r0o(g)— = -2zcos(rago) 

In other words 

n dq sininq) 
^'^ / ^ x{q)-l{q') = = '^(sin(ng)) 



Therefore 



In X space. 



-V2 - y) 

It follows that the eigenvalue equation 



^ sin(ng(y)) 

ay— — = — cos(?2g(x)) 



d /"^ dy7C(f)o{y)fn{y) _ , 



dx i_V2 (X - 

has solutions 



sin(ng(x)) sin(ng(x)) AT ^ n 

fn{x) = 7 = , x{q) = -V2cos(g), e„ = n 

vr0o V2sin(g(x)) 

This follows from the observation that in terms of time of flight coordinates, 
the above spectrum equation takes the form 

f dq' 7T(j)o{q')fn{q') ., , . w . ^^ 

J - x{q) - x{q') ^-^\9.\i-M9)fnW)) 

Concerning the first term in ()57|) . we have already seen for the main text 
that it can be obtained straightforwardly from the result (eqn. (|iup) 

dz-^ = X. (59) 
[x — z) 
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This equation is solved by the well known methods of ref. ^Hl- We point 
out that first term of (|57|) can also be obtained in general by considering the 
integral 



along the contour described above. There is now a contribution from " infin- 
ity", and one obtains the result that the above integral equals —1 



7.4 Appendix D: Three impurities 

For three impurities, we have 

^(fc,3) = [ dk2 r dkiTr{Be'^'^ Be'^^'-^^'^^ Be'^^-^^^^'^) 
Jo Jo 



(60) 



-3J:,,,4V^BV),,iV+BV),kiV+BV) 



ki 



{\j — — \k) 



and 



V^(x,3) = -3J:^,,,kiV-^BVy,,iV+BV),kiV+BV)t. 



After some algebra, one obtains 



6{x — Xj) 



{x - Xi){x - Afc) 



(61) 



n{ko,0 : fc,3) 



-ko^,,,4V+BV),,{V+BV),k{V+BV)k^ 
1 



(Aj — Aj)(Aj — Afc) 



\{Xi - Afc)2(A, - A,) ^ (A, - A,)2(A, - A,) 
^ (A, - A,)2(A, - A,) (Afe - A,)2(A, - A,) 



and 
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(l{x,0 :y,3) = j:ij,t{V*BV),i{V+BV)jt(V+BV)t, X 

+ ^^-^ (A, - - A.) ^ + ^^-^ (A. - - A,) I 
The f2(a;, : y,3) term in ()13p takes the form 

. S{y-\i) din J d 



- 5(a;-Ai)5(?/- Ai)(- — — — —)d^ 



{y - A,)2 {y - \) {y - {y - \^)' d^{x, 0) dil^^y, 3) 
f 5{x - Xj)5{y - Xj) ^ 5{x ~ Xk)6iy ~ Xj) ^^ d\nJ d 



iX,-yny-Xk) (Xk-ynv-X,) > 0) 3) 

For the harmonic oscillator potential, we use the results (fTT|) and ()10p , so 
that we can write 



'dilj{x,0) J {x — z) 

Then 
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1 /",/", , dlnJ d 

ax / dyil[x, : y,3) 



2 J y ^ V ' • ^' ' dilj{x,0)dij{y,3) 

6{x-Xi) f , Mz) r, d 



+ 3 j dyJ:,,,4V+BV),,{V+BV),k{V+BV)k^6{y - X,] 



d 



dip{y,3) 



/ X{ 1 Aj 1 Xj Afc 

^[y-X^Y{y-X,) + {y - Xjy {y - A,) " (A, - y^y - A^) ~ [X^-yfiy-X,) 

dx f dz f'^^ 3)(9^— — r~2 /" 3)(9i: 



(x — 2;) ' dip{x, 3) 7 ' dil){x, 3) 



Again, the first term above cancels the similar term in |T3jl , and we obtain 
for the quadratic hamiltonian in the 3 impurity sector: 

- j dxjdz ^^^^^ 2 J dxiJ[x, 6) ^^^^^ 

(62) 

This is again a shifted Marchesini-Onofri operator. The spectrum and eigen- 
functions of this operator are 

, ,_o sm(nq) 

Wn = n-3 ■ (t>n = . , \ ^=1,2,... 

V2 sm(gj 

Adding the contribution from the Tr{Bd/dB) term of the Hamiltonian we 
obtain 

, ,_o sinfno) , , 

Wn = n ; 0r'= ^ \ \ ; n = l,2,... (63) 

V 2 sm(gj 
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